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ABSTRACT 

We present an expression for the nonlinear evolution of the cosmological power 
spectrum based on following Lagrangian trajectories. This is simplified using the 
Zel'dovich approximation to trace particle displacements, assuming Gaussian initial 
conditions. The model is found to exhibit the transfer of power from large to small 
scales expected in self-gravitating fields. Some exact solutions are found for power-law 
initial spectra. 

We have extended this analysis into redshift-space and found a solution for the 
nonlinear, anisotropic redshift-space power spectrum in the limit of plane-parallel 
redshift distortions. The quadrupole-to-monopole ratio is calculated for the case of 
power-law initial spectra. We find that the shape of this ratio depends on the shape 
of the initial spectrum, but when scaled to linear theory depends only weakly on the 
redshift-space distortion parameter, (3. The point of zero-crossing of the quadrupole, 
fco, is found to obey a simple scaling relation and we calculate this scale in the 
Zel'dovich approximation. 

This model is found to be in good agreement with a series of TV-body simulations 
on scales down to the zero-crossing of the quadrupole, although the wavenumber at 
zero-crossing is underestimated. These results are applied to the quadrupole-monopole 
ratio found in the merged QDOT plus 1.2 Jy IRAS redshift survey. Using a likelihood 
technique we have estimated that the distortion parameter is constrained to be p > 0.5 
at the 95% level. The local primordial spectral slope is not well constrained, but the 
likelihood analysis suggests n « — 2 in the translinear regime. The zero-crossing scale 
of the quadrupole is ko = 0.5 ± 0.1/iMpc -1 and from this we infer the amplitude of 
clustering is as = 0.7 ± 0.05. 

We suggest that the success of this model is due to nonlinear redshift-space effects 
arising from infall onto caustics and is not dominated by virialised cluster cores. The 
latter should start to dominate on scales below the zero-crossing of the quadrupole, 
where our model breaks down. 

Key words: Cosmology: theory - large-scale structure of Universe 



1 INTRODUCTION 

The Newtonian analysis of linear growth of perturbations in 
an expanding universe is a well understood problem (Pee- 
bles 1980; Efstathiou 1990). The extension of this into the 
nonlinear regime has proven more difficult due to the strong 
mode coupling that arises in gravitational collapse, and most 
progress has been made through the use of iV-body simula- 
tions. 

Actual observations of galaxies in redshift space are fur- 
ther complicated (and made more interesting) by redshift 



distortions, caused by peculiar velocities of galaxies along 
the line-of-sight. Again, the linear problem is relatively well 
understood. Kaiser (1987) has shown that linear redshift dis- 
tortions take their simplest form when expressed in Fourier 
space, at least if structure is far from the observer so that the 
distortions are essentially plane-parallel. Here the redshift- 
space Fourier modes are related to the real-space ones by 

8 s (k) = (l + Pfi 2 k )5(k) (1) 

where a superscript s denotes a redshift-space quantity, [ih 
is the cosine of the angle between the wavevector k and the 
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line of sight, and f3, the redshift distortion parameter, is the 
dimensionless growth rate of growing modes in linear theory, 
which is related to the cosmological density £1 by (Peebles 
1980) 



(2) 



in the standard pressureless Friedmann cosmology with 
mass-to-light bias b. It is through measuring the distortion 
parameter j3 that one hopes to measure the cosmological 
density parameter, fi. 

Kaiser's formula (Q) is valid only in the linear limit, 
and for plane-parallel distortions, neither of which approxi- 
mations is well satisfied in reality. Several authors have now 
addressed the problem of generalising Kaiser's formula to 
the case of radial distortions, while retaining the assumption 
of linearity (Fisher, Scharf & Lahav 1994; Heavens & Taylor 
1995; Fisher et al. 1995; Zaroubi & Hoffman 1994; Ballinger, 
Heavens & Taylor 1995; Hamilton & Culhane 1996). 

The main purpose of the present paper is to extend the 
analysis of redshift distortions into the nonlinear regime, re- 
taining the plane-parallel approximation for simplicity. Our 
approach to the problem is motivated by the consideration 
that the density in redshift space may appear highly nonlin- 
ear even when the density in real space is only mildly non- 
linear. For example, a region which in real space is just turn- 
ing around, a mildly nonlinear condition, appears in redshift 
space as a caustic, a surface of infinite density, which is thor- 
oughly nonlinear. This leads us firstly to work in Lagrangian 
space (Section ^]), and secondl y to adopt the Zel'dovich 



2j & |). The Zel'dovich 



(1970) approximation (Sections 

approximation is in effect linear theory expressed in La- 
grangian space, inasmuch as it approximates the trajectories 
of particles as straight lines with (comoving) displacements 
growing according to linear theory. 

Our approach follows that of Taylor (1993), who studied 
the nonlinear evolution of the power spectrum. Comparable 
approaches have been used by Bond & Couchman (1987) 
to evolve the galaxy angular correlation function, by Mann, 
Heavens & Peacock (1993) to evolve the real-space corre- 
lation function of clusters, and by Schneider & Bartelmann 
(1995) to evolve the real-space (unredshifted) power spec- 
trum. Our approach differs from that of Hivon et al. (1995), 
who applied a perturbation expansion in Lagrangian space 
to second order to calculate the redshift-space skewness. 

For simplicity, we assume throughout this paper that 
the density field is unbiased, 6 = 1. Generally, the effect of 
evolution, if continuity is assumed (which is a fundamental 
assumption of this paper) , is to tend to drive the bias factor 
towards unity. Continuity implies that the ratio of galaxy 
to matter density (1 + 5)/(l + 8m) remains constant in La- 
grangian elements, and if (somehow) a linear bias 8 = &<5m 
is established at some early time when 5 and 5m are both 
small, then the ratio of galaxy to matter density must be 
close to unity, (1 + 6<5m)/(1 + <5m) ~ 1- It follows that the 
bias will be close to unity at later, nonlinear epochs when Sm 
is no longer small. Conversely, if bias is in fact important at 
the present, nonlinear epoch (as may well be the case), then 
it must be that the assumption of continuity must break 
down in the not too distant past. Indeed, it may be that 
continuity is violated on an ongoing basis. In the present 



paper we choose to ignore this thorny problem, and simply 
assume an unbiased density field. 

We begin in Section H by deriving equations which re- 
late the Lagrangian and Eulerian descriptions of density, 
and we calculate a general expression for the evolution of 
the power spectrum. We invoke the Zel'dovich approxima- 
tion and Gaussian initial conditions and consider some of 
the general features of the resulting nonlinear power spec- 
trum. Some exact solutions for initially power-law spec- 
tra are derived. In Section ^ we calculate the power spec- 
trum in redshift-space, again in the Zel'dovich approxima- 
tion, and obtain a number of analytic and numerical results 
for the observationally interesting ratio of quadrupole-to- 
monopole power. In Section ^, we compare the predictions 
of the Zel'dovich approximation with JV-body simulations, 
and we apply our findings to analyse the quadrupole distor- 
tion measured in the QDOT plus 1.2 Jy redshift survey. We 
summarise our conclusions in Section ra. 



2 REAL-SPACE POWER SPECTRA 

2.1 Clustering in the Lagrangian Frame 

In the Lagrangian theory of fluid mechanics, the central vari- 
able is the displacement vector field, rather than the density 
field of Eulerian space. This displacement field, £(q, t), maps 
particles from their initial Lagrangian coordinate, q, to the 
Eulerian-space coordinate, x, via the relatiorj^] 

x(q,t) = q + Z(q,t) (3) 

where £(q,t) is the integral of the velocity field along the 
world-line of the particle: 



£(q,t) = J dt' v[x(q,t'),t'] 



(4) 



The advantage of the Lagrangian formalism for gravita- 
tional collapse was first pointed out in the seminal work 
by Zel'dovich (1970). Because the coordinate system of 
Lagrangian space is nonlinear, moving with the particles 
themselves, this approach takes the analysis of perturbation 
growth of the density field into the nonlinear regime while 
conserving mass-density. 

However, the most useful observable quantity in cosmol- 
ogy is the Eulerian density field, p(x), inferred either from 
variations in the galaxy distribution, or indirectly from dis- 
tortions in the Hubble flow, gravitational lensing, or from the 
microwave background. Continuity, along with the assump- 
tion of uniform initial density, implies that the relation be- 
tween the overdensity 5(x) = [p(x) — p]/p in Eulerian space 
and the Lagrangian displacement field $,{q) is 



5{x) 



d qS D [x q £(q) 



(•») 



where 5d(x) is the Dirac delta function. The Fourier trans- 
form 



S(k) 



d 3 x e lk x S(x) 



(6) 



* Throughout this paper we take the coordinates x and q to 
be comoving coordinates, defined relative to the general Hubble 
expansion of the Universe. 
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of the overdensity S(x), equation (|^), is 
S(k) = / d 3 qe ikq (e ik( - 1) , 



(7) 



which relates the Fourier space Eulerian density field 5 (As) 
to the Lagrangian displacement field £(q). 

The correlation function of Fourier modes defines the 
power spectrum P(k): 

(5(k)5*(k')) = (2nfP(k)8 D {k - k') , (8) 

where the angular brackets denote ensemble averaging and 
the Dirac delta-function arises from translational invariance. 
Transforming to differential and centre of mass coordinates, 
we find that as a result of translational invariance only the 
differential terms survive. We can then express the power 
spectrum by the integral equation 



P(k) = / dV 9 ((e 6 ) - 1) 



(9) 



where the term (e lfe,A ^) is identified as the generating func- 
tion of the differential displacement vector field, A£(q), of 
points separated by distance q = q 1 — q 2 : 



A£(<?)=£(<h)-£(q 2 



(10) 



Equation (|9j) for the power spectrum is valid at all times, 
not merely in the linear regime. However, in order to solve 
this completely one needs to have an expression for the dis- 
placement generating function and its evolution. In practice 
we shall have to approximate this. In the next Section we 
use equation (^) to calculate the evolution of the power spec- 
trum in the Zel'dovich approximation, under the assumption 
of random Gaussian initial conditions. 



2.2 Evolution of the power spectrum in the 
Zel'dovich approximation 

A careful analysis of the Zel'dovich approximation has been 
given by Hui & Bertschinger (1996), who consider it in the 
context of local theories of gravity. For our purposes we only 
need the results that the particle displacement field scales 
according to linear theory: 



(2tt) 3 



d 3 fc Aj L ( fc ) e - 



ik.q 



(11) 



where £l(<z) is the linear displacement field, and Sh(k) the 
linear density field, defined at some suitably early time. 
This approximation is then a first order theory in the dis- 
placement vector field, extrapolated to arbitrary later times. 
By construction the Lagrangian coordinate system preserves 
mass, but in the Zel'dovich approximation at the expense of 
not satisfying the Euler equation. In the Lagrangian frame 
the density field is usually evolved along fluid streamlines, 
5[q(x), t], but equation (El) allows one to calculate the den- 
sity in Eulerian coordinates. 

The assumption of a Gaussian initial fluctuation field 
implies that the initial displacement field, £ L , was Gaussian. 
The Zel'dovich approximation ( |ll| ) further implies that the 
displacement field £(q) remains Gaussian for all time. Equa- 
tion (^) for the power spectrum involves the Lagrangian gen- 
crating function of the differential displacement field, which 



for a Gaussian vector field is 

(e ife - A «> = exp(-fc l fc J (A^AC 3 )/2) 

= exp(-kikj [ipij (0) - Vij (q)]) 



(12) 



where ipij(q) = (£i(<7i)£j(<72)) i s the displacement corre- 
lation function. In the Zel'dovich approximation, the dis- 
placement correlation function ipij grows as the square of 
the growth factor D(t) in the nonlinear as well as linear 
regimes: 

ipii(q,t) = [D(t)] 2 ^ ij (q) . (13) 

The linear correlation function iphii(q) of displacements is 
related to the linear power spectrum Ph(k) by 

ipu 3 {q) = (£ L i(<?i)£ Lj (<? 2 )) 

= (2^/ d3fc ^ PL(fe)e " ife - 9 
which at zero separation, q = 0, is 



(14) 



^y(O) = {HSU = 



3(2tt) 3 



d J kk~ 2 P L (k) 



(15) 



where (£ 2 ) is the one-dimensional dispersion of displace- 
ments. 

Inserting the Gaussian generating function ( |l2| ) into 
equation (Bj) for the power spectrum yields (Taylor 1993) 



P(k) = / d 3 qe ik - q {e- k ' k ^^ (0) -^ (q)] -1} . 



(16) 



This equation ( Jl6| ) , along with equations ( |l3| and ( |l4| ) , pro- 
vides a nonlinear mapping between the initial and evolved 
power spectra, where the degree of evolution is controlled 
by the linear amplitude, D(t), of perturbations. 

The unit term on the right of the integrand of equa- 
tion @, which arises from subtraction of the mean den- 
sity, yields a delta-function at k — on integration. The 
unit term can be ignored for k ^ 0, and we tacitly drop it 
hereafter. Nevertheless, it can help in numerical integrations 
to retain the unit term when k is small (linear regime). 

So far we have not assumed that the density field is 
isotropic, although we have assumed it is homogeneous. In 
Section g we will apply the power spectrum equation (|l^) in 
redshift space, where the density field is not isotropic. In the 
remainder of the present Section we assume that the density 
field is unredshifted and isotropic. 

For an isotropic density field, the displacement corre- 
lation function ipij(q) resolves into irreducible components 
parallel and perpendicular to the pair separation q: 



where 



(?) = 
and 

i>±(q) 



[D{t)f 
2tt 2 



[D{t)f 
2tt 2 



dfeP L (fe) 



dfcP L (fc) 



jo(kq) - 2 



ji(kq) 



ji (kq) 



kq 



kq 



(17) 



(18) 



(19) 



For a general potential flow, as here, the parallel vari- 
ance tpii is related to the perpendicular variance tj}± by 
(Monin & Yaglom 1975; Gorski 1988) 



© 0000 RAS, MNRAS 000, 000-000 



4 A. N. Taylor & A. J. S. Hamilton 



^11(9) = 



d#±(g) 

dq 



(20) 



which can also be derived from equation ( p^ ) with the def- 
initions (Jl7|) . It is convenient to define the differential dis- 
placement covariances ip+ and xp- by 



^+(q) = ^_l(0) - tl>±(q) 
%j)-{q) = ip x (q} - -i/>||(g) 
which are related by 
qdip+(q) 



Aq 



(21) 
(22) 

(23) 



The (isotropic) power spectrum ( |l6| ) can then be written in 
terms of 4i+ and ip- 



P(k) = d\e 



\kqy,-k (j/i + + !/>_/J ) 



(24) 



where /x = k.q. Integrating equation ([mJ) over the azimuthal 
angle of q about the wavevector k is trivial, and further 
integration over fi yields 



>F(kq,k 2 ip- 



P{k) = Tie \ 1-Kq £ dqe 
Jo 

where F(A, B) is the function 
_i/2 / 42, 



1/2 



_L4_ 

2SV2 



(25) 



(26) 



with erf(z) the error function (Abramowitz & Stegun 1968) 



erf (2) 



2 

^72" 



dte" 



(27) 



The error function erf(z) for complex z is available in pro- 
grams such as Mathematica. Equation (JSHJ) reduces the eval- 
uation of the power spectrum to a single integral ove r sep- 
aration q. Excepting some simple cases (Section 2.4), this 
integral must be done numerically. 

The integrand in equation (BO) oscillates rapidly when 
kq is large and real, which can make the integral difficult 
to evaluate (Schneider & Bartelmann 1995). One way to 
resolve the difficulty is to continue the ip±(q) analytically 
into the complex plane, and to shift the path of integration 
over q into the upper complex plane, which converts the 
oscillations into exponential decay. Some experimentation 
may be required to determine a good integration path, along 
which the integrand is well-behaved. For example, in the case 
of an initially Gaussian power spectrum, Section 2.3 below, 



we found it satisfactory to integrate first along the real axis 
to q ~ min(2D, 1/k), and then to complete the integration 
along a straight line slanted upwards in the complex plane, 
at angle 7r/6 from the real line. Variants on this path, with 
or without the initial segment along the real axis, and with 
other slant angles, worked in other cases. 



2.3 Example: a Gaussian power spectrum 

As an example, in Figure [I] we show the Zel'dovich evolution 
of a power spectrum which is initially gaussian, 



O0.1 




100 



Figure 1. The evolution of the dimensionless amplitude Al(&) = 
fc 3 / 2 exp — fc 2 /2 in the Zel'dovich approximation for random Gaus- 
sian initial conditions. The solid line is the nonlinear power spec- 
trum and the dotted line is the linear initial spectrum. 



The quantity plotted is the dimensionless amplitude A(fc), 
the square root of the dimensionless power spectrum defined 
by (Peacock 1992) 



A 2 (fc) 



k s P{k) 
2tt 2 



(29) 



Ph(k) — exp(— k 



(28) 



Figure |l] shows clearly the two main features of nonlinear 
growth, the amplification of small scale power from the for- 
mation of caustics, and the restriction on the growth of in- 
termediate scale power imposed by the positive density con- 
straint in voids. 

Figure |l| also illustrates the limitations of the Zel'dovich 
approximation at later stages. While in reality the small 
scale density contrast increases following the collapse of 
structures, in the Zel'dovich approximation structures ex- 
pand back outwards after collapse. Thus in Figure ^] the 
amplitude A(fc) reaches a maximum of order unity as a func- 
tion of time, and thereafter declines. In reality, the ampli- 
tude A(fc) would be expected to increase monotonically with 
time. 



2.4 Solutions for initial power— law spectra 

For initially power-law spectra, Ph(k) oc k n , convergent re- 
sults are attained in the Zel'dovich approximation for power 
spectrum indices in the range —3 < n < — 1. For n < —3, the 
differential acceleration between particle pairs, hence the dif- 
ferential displacement field, receives divergent contributions 
from large scales, while for n > — 1, the acceleration diverges 
on small scales. In the latter case, nonlinear processes are 
expected to intervene so as to truncate the divergence, but 
we do not consider such truncations here. 

For initially power-law spectra the differential displace- 
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ment covariances are power laws; ip± oc q n . It is conve- 
nient here to normalize the growth factor, D(t), so that 

iP + {q)=D 2 q~ n - 1 (30) 

(a normalization based, for example, on the power spectrum 
is less convenient because it leads to divergence in the dis- 
placement covariance for n = — 3 and —1). The relation ( |23| ) 
implies 



iMg) = (-n-l)V+(g) 



(31) 



With the normalization (|30j), the initial linear power spec- 
trum Ph(k) is (note Ph(k) is defined without the factor D 2 , 
so -P(fc) = D 2 P\ J {k) in the linear regime) 



PL(fe) = 4tt(2-ti) sin[(-n-l)7r/2]r(l-n)fc ri 



(32) 



which is positive over the interval —3 < n < —1, going 
through zero at n = —3 and —1. 

In the Zel'dovich approximation, initially power law 
power spectra evolve in a self-similar fashion: 



A 2 (k,D) = A 2 (D 2 k 



2,n + 3x 



(33) 



The self-similar evolution of scale-free initial conditions in 
the case of a flat, Q = 1, universe is familiar from the 
BBGKY hierarchy (Peebles 1980). In the Zel'dovich approx- 
imation, the property of self-similarity extends also to non- 
critical Universes. 

We now give exact solutions for the cases n = —1, —2, 
and —3 in the Zel'dovich approximation. 



(i) n - 



In the limiting case n — > — 1, the displacement field has an 
incoherent Gaussian nature. This arises because the incipi- 
ent divergence of acceleration on small scales causes particles 
at finite separations to acquire essentially random velocities 
relative to each other. But if particles in an initially uniform 
field move about randomly, the resulting field continues to 
remain uniform. Therefore, the power spectrum should re- 
main infinitesimal even when the displacement covariance 
grows to a finite value. Mathematically, the linear power 
spectrum (^) goes to zero in the limit n — > — 1. 

In reality, nonlinear evolution will in this case invali- 
date the linear extrapolation of displacements assumed by 
the Zel'dovich approximation. Nevertheless, it is interesting 
to consider n = — 1 as a limiting case. Defining the small 
positive quantity v by 



v = — n — 1 



(34) 



which tends to zero as n — > — 1, one finds that for ip+ = 
D 2 q v , equation (jsol), the dimensionless power spectrum, A 2 , 
evolves as 



A 2 (k,D) = 3vD 2 k 2 e~ 



(35) 



The power spectrum (j35|) is that of an n = -1 Gaussian 
field, modulated by a random Gaussian distribution of dis- 
placements with a one-dimensional one-point dispersion of 
D. 

Although the power spectrum is itself infinitesimal, 
the ratio P| / Pq of quadrupole-to-monopole redshift power, 
which we shall consider in more detail later, is well-defined 
and finite in the limit n — > — 1. 



(ii) n = -2 

For an initial power-law spectrum of index n = —2, the 
integral ( |24| ) can be done analytically. Here ip+ = tp- = D 2 q, 
and the dimensionless power spectrum evolves as 



A 2 (k,D) = 



lQD 2 k 



tt(1 + 4D 4 fc 2 ) 2 



1 + 



37rZ> 2 fc 



(1 + 4P> 4 fc 2 ) 1 / 2 



(36) 



In the linear regime, A 2 = 16D 2 k/ir, while in the highly 
nonlinear regime A 2 oc (D 2 k)~ 3 . 



(iii) n - 



In the limiting case n — > —3, differential accelerations be- 
tween pairs are dominated by tidal contributions generated 
on large scales, inducing motions which appear locally like 
anisotropic Hubble flows. The displacement covariances are 
i/>+ = D 2 q 2 , ip- — 2ip+. The integral (^) can again be done 
analytically, yielding 



A 2 (k,D) 



e -l/{12D 2 ) 

(12n) 1 / 2 D 3 



(37) 



Here the power spectrum remains a power-law, P(k) oc fc~ 3 , 
at all times. As argued by Schneider & Bartelmann (1995), 
the persistence of the k~ 3 form results from the dominance 
of caustics in the density distribution. The amplitude ini- 
tially increases, from an exponentially tiny value, as more 
regions reach the point of collapse, but then reaches a max- 
imum, at D 2 = 1/18, and subsequently declines as D~ 3 
as regions pass through the point of collapse and continue 
streaming outward. 

We note that equation (^) also gives the asymptotic 
form of the evolution of the power spectrum in the Zel'dovich 
approximation on small scales whenever the spectral index 
is less than —3 on small scales. This is true for example in 
the case of the Ga ussi an power spectrum Phik) = e~ k con- 



sidered in Section 2.3. In this case, rms differential displace- 
ments between pairs are proportional to their separation on 
small scales. Define 



lim 

9^0 



D 

60tt 2 



dkk 2 P h (k) 



(38) 



Then for large k (larger than the wavenumber at which the 
spectral index passes through n = —3), 



A (k,D) 



-l/(12a 2 ) 



(39) 



(127r)VV 3 

which is the same as equation (pS?]), with a in place of D 



3 REDSHIFT— SPACE DISTORTIONS 

In Section || we considered some of the properties of the den- 
sity field evolving in the Zel'dovich approximation. Here we 
turn our attention to the properties of such a density field 
when viewed in redshift -space. We assume the Zel'dovich 
approximation throughout this Section; in Section 5 we will 
test the validity of the approximation with iV-body simula- 
tions. 
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3.1 Zel'dovich approximation 

In redshift -space, the observed radial position of a galaxy 
is given by its radial velocity, which measures not only the 
uniform Hubble flow, but also the peculiar motion of the 
galaxy. Thus the position s of a galaxy in redshift space 
appears shifted from its true position x by its radial peculiar 
velocity^ 

s = x + [x.v(x)]x . (40) 

In the Zel'dovich approximation, the peculiar velocity v is 
related to the displacement £ at all times by the linear re- 
lation v[x(q),t] = £(q,t) — f3£(q,t), where f3, equation (|^), 
is the dimensionless linear growth rate. As discussed in the 
introduction, we assume an unbiased density field, 6 = 1. 
The position s of a galaxy in redshift space then appears 
displaced from its initial Lagrangian position q (this initial 
position q is the same in both real and redshift space) 



s = q + e(q) 



(41) 



by a redshift displacement £ s (q) which in the Zel'dovich 
approximation is 

V(q) = t(q)+P[&4(q)]&- (42) 

Equation ( ffl| ) shows that the displacement £ s of a galaxy 
in redshift space is stretched by a factor (1 + /3) along the 
radial direction. The overdensity 8 3 in redshift space is then 
(compare equation M) 



d q5 D [s q £ s (q) 



(43) 



The radial character of the redshift distortion /3(x.£)x 
complicates taking the Fourier transform of the overdensity. 
However, matters simplify if the structure being observed is 
far away from the observer, so that the distortions are ef- 
fectively plane-parallel. In the plane-parallel approximation, 
the redshift displacement £ s becomes 



(44) 



where the line-of-sight direction z is taken to be fixed. The 
Fourier transform of the redshift density, equation ([l3|), is 
then 



s s (fe) 



J d 3 qe ik q (e ik iS 



(45) 



It is convenient to introduce the vector K, which is the 
wavevector k stretched by (1 + /3) along the z-axis: 



K = k + (3(z.k)z 

whose magnitude is 

K = k(l + 2/3/4 +/3 2 Mfe) 1/2 



(46) 



(47) 



with fik the cosine of the angle between the wavevector and 
the line of sight 

Mfe = z.k . (48) 
In terms of the vector K, equation (f45J) can be rewritten 



t From hereon we work in distance units of the velocity, where 
the Hubble parameter Hg = 1. 



5 s (k) = / dV V - 1 



(49) 



Just as the displacement £ s of a galaxy in redshift space 
is stretched by a factor (1 + /3) along the line-of-sight axis z, 
equation (Q), so also is the displacement correlation func- 
tion ipij = (£f£^) in redshift space stretched by (1 + 0) 
factors along each z-component of the tensor: 



ipij = {5 ik + f3ziZ k )ipki(5ji + P&jZi) 



(50) 



If the unredshifted displacement field £ is Gaussian (as is 
true for an initially Gaussian field evolved in the Zel'dovich 
approximation, Section 2.2), then the redshifted displace- 
ment field £ s is also Gaussian. 

The derivation of expression (^) for the evolved power 
spectrum in Section 2.2 did not assume isotropy, and re- 



mains valid in redshift space. Thus the redshift power spec- 
trum P s (k) is given by equation (^6|) with the displacement 
correlation function taken to be the redshift space quantity 



P s (k) 



a qe \e * 3 1 ^ K ' »j ■ " — 1} 



(51) 



As in equation (^), the unit term in the integrand can be 
dropped for k ^ 0, although again it helps in numerical 
integrations to retain the term when k is small. Symmetry 
about the line of sight z implies that the redshift power 
P s {k, Hk) is a function of the magnitude k of the wavevector 
k and of the cosine \Xk of the angle between the wavevector 
and the line of sight. 

The relation ( jHoj ) between the redshifted and unred- 
shifted displacement correlation functions ipfj and ipij allows 
the combination kikji^lj in the integrand of equation (jHTJ) 
to be written conveniently in terms of the vector K , equa- 
tion (^), as kikjipij = KiKjipij. Then equation dEl] ) for the 
redshift power spectrum becomes 



P s (fc) = 



(52) 



where now n = K .q, and tl>±(q) are the (unredshifted) quan- 
tities given by equations (Gil) and (E2|). Equation ( p2[ ) is 
the redshift space counterpart of equation (E4). Integrating 
equation (^) over the azimuthal angle of q about the vector 
K yields a Bessel function Jo 

P s (k)= ( 2nq 2 dqdfie [kq ^- I<2 ^+ + ^~» 2) J (kqs v / T^)(53) 



where s and c are 



,1/2 



c=(l 



(1 + 20»l + pfi 
s 2 ) 1/2 = 



2..2 V /2 



(1 + 2^1+/?VD 1/2 



(54) 



(55) 



Excep t fo r the interesting case s = 
in Section 



3.4 



pursued further 
cannot be reduced further to 



equation 

any simple analytic form. Moreover, the oscillatory character 
of the integrand makes it difficult to evaluate numerically. 
Further details on the evaluation of the integral (^) are 
given in Appendix A. In Section 3.5 we adopt a different 



approach to evaluating the redshift power spectrum, which 
is to smooth it. 
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3.2 Redshift power spectra for initial power— laws 

We now give analytic solutions for the evolution of the red- 
shift power spectrum in the Zel'dovich approximation for 
initially power-law spectra Pl oc k", for the cases n = — 1 
and n = —3. 

In Section ^4| we saw that the n = — 1 spectrum gave 
rise to an incoherent random Gaussian velocity dispersion, 
resulting in a simple damping term for the nonlinear power 
in the Zel'dovich approximation. Extending this to the red- 
shift domain, we can see that the effect will be an additional 
damping term in the line of sight direction — an approxima- 
tion previously suggested by Peacock (1992) to model the ef- 
fects of virialised clusters. With the normalisation discussed 



in Section 2.4, we find that the redshift power spectrum is 
A s2 (k,D) = 3vD 2 k 2 (l+(3/il) 2 e~ D2K2 . (56) 

This is the same as the Peacock smoothing with a 2 — 
D 2 (20 + 2 ), where a v is the one-dimensional one-point 
velocity dispersion. Notice that we also recover the linear 
Kaiser boost factor, (1 + /3/x 2 .) 2 . Hence we conclude that 
this form has some justification beyond its phenomenologi- 
cal foundations for scales where the spectral index is ~ —1, 
and arises as a special case in the model presented here. 

For n — —2, there appears to be no analytic solution for 
the redshift power spectrum, unlike the unredshifted case, 
equation (p6|). 

In the case n = — 3, we find the analytic solution 

A a2 (k,D) ' 



(127r)V2D3(l + 2pnl + /3 2 /4) 3 /2 



1 + 2/3 M 2 + 3/3V1 - 2/3 2 /4 



exp 



12.D 2 (l + 2/V fc +/3 2 /4) 2 



(57) 



which is the redshift counterpart of the unredshifted solution 
©• 

3.3 Multipole analysis of the redshift power 
spectrum 

The effect of redshift distortion on the statistical properties 
of the density field can be best interpreted by a multipole 
expansion of the redshifted power spectrum. The anisotropic 
spectrum can be expanded in Legendre polynomials 



p s {k) = ^pi(k)vt(,iik 



(58) 



where Pe(fik) is the i Legendre polynomial and fik is again 
the cosine of the angle between the wavevector k and the 
line of sight. The multipole moments of the power are given 
by 



Pe{k) = 



2£ + l 



dp fe P s (k)Ve{n k ). 



(59) 



Hamilton (1992) suggested using the ratios of the mul- 
tipoles as a way of measuring 0, expanding the correlation 
function from Kaiser's linear theory. Cole, Fisher & Wein- 
berg (1994) developed this to study the anisotropy of the 
power spectrum, where nonlinear effects could be avoided 
by limiting the analysis to small wavenumbers. The distor- 
tion parameter, 0, could then be estimated from the ratio 
of monopole to quadrupole power: 



Pg{k) 



IP + IP 2 
l + |/3 + i/3 2 



(60) 



In the linear regime each observable mode contributes to an 
independent estimate of 0. 

Nonlinear effects modify the quadrupole to monopole 
ratio from the linear value (^). In the following Section 
we derive an approximate expression for this ratio in the 
Zel'dovich approximation, and in subsequent Sections we 
present numerical results for the ratio in the Zel'dovich ap- 
proximation, in iV-body simulations, and in observations. 

3.4 An approximation to the quadrupole to 
monopole ratio 

In this Section we derive an approximation, equation (^S4j), 
which relates the ratio P| / Pq of the quadrupole to monopole 
redshift power spectrum to the evolution of the unredshifted 
power spectrum in the Zel'dovich approximation. 

In equation (j^j for the redshift power spectrum, notice 
that the Bessel function Jo{kqs^Jl — fi 2 ) in the integrand 
equals one whenever s is zero, which happens in the limit 
— > (equation [Q) and also in the cases fj,u = or 1 
where the wavevector is perpendicular or parallel to the line 
of sight. Moreover in the nonlinear limit where k — > oo, the 
K 2 terms in the exponential cause the integrand of equation 
( p3| ) to decay rapidly as q increases and again the Bessel 
function is sensibly equal to one (as is also the factor e lkqcti ) 
wherever the integrand is non-negligible. In all these cases, 
the redshift power, equation (^), reduces to the same form 
as its unredshifted counterpart, equation (pi|), but with ip± 
in the latter replaced by (K/k) 2 ip±, where K/k = 1 + 0a 2 , 
(equation ^t|) given that s = and c — 1 (equation fig]). 
In other words, the redshift power equals the unredshifted 
power evolved forward by l + /3pi| in the growth factor D(t): 



P 8 (k,D)^P[k,(l+0^ 2 k )D] 



(61) 



As it happens, this same equation (^SlJ) is also valid in 
the linear regime, according to the usual linear relations 
P(k,D) = D 2 P L (k) and P s (k,D) = (1 + 0^ 2 k ) 2 P(k, D). 
Thus approximation (^l|) is valid (a) in the limit — > 0, or 
(b) if fik = 0, i.e. if the wavevector k is perpendicular to the 
line of sight, or (c) if /i*, = 1, i.e. if the wavevector k is paral- 
lel to the line of sight, or (d) in the highly nonlinear regime, 
or (e) in the linear regime. Evidently the approximation is 
of some generality. 

Let us consider more closely the limiting case /3 0, 
for which the approximation ( |6l"| ) becomes exact (within the 
context of the Zel'dovich approximation). Expanding equa- 
tion (|l|) as a Taylor series for small gives 



P s (k,D)=P(k,D)+0fi 2 k D 



dP(k,D) 
dD 



09 -0) 



(62) 



Here the redshift power spectrum is the sum of a con- 
stant term and a term proportional to n%, which implies 
that the redshift power spectrum is a sum of monopole and 
quadrupole terms, for small 0. The ratio of quadrupole to 
monopole power for small is, from equation (p2J), 



40 d\nP(k,D) 
~3~ 



03-0) 



P s (fc) 3 <91n(D 2 ) 
Amongst other things, equation (r 



(63) 

predicts that the 
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quadmpole redshift power P£ (k) goes through zero where 
the unredshifted power P[k,D(t)] reaches its maximum as 
a function of time in the Zel'dovich approximation. 

We now propose a generalisation of the result ( |63| ) to 
arbitrary 0, which is consistent with and to some extent 
motivated by approximation (|6l|), and which accords with 
the numerical results of Section j3.6| below, which are further 
supported by the iV-body simulations in Section [O]. The nu- 
merical results indicate that, at least in the case of initially 
power law power spectra, the quadrupole to monopole ratio 
Pi I Pq appears to satisfy a simple scaling law with 0, such 
the shape of the ratio is insensitive to 0, while its amplitude 
is proportional to the usual linear ratio (strictly, the numer- 
ical results are for smoothed power spectra, but we assume 
that the results hold also without smoothing). We further 
find that the zero-crossing of the quadrupole scales with 
such that D 2 ~ 1/(1 + 0). These empirical results, combined 
with equation 



suggest the general approximation 
f/3+f/3 2 dlnP[k,(l+0)^ 2 D] 



(64) 



Pg(k) " l + ain^) 

The approximation ( pi] ) relates the quadrupole to monopole 
ratio P| /Pq of the redshifted power spectrum to the evolu- 
tion dlnP/dlriD of the unredshifted power spectrum in 
the Zel'dovich approximation. Although we have numerical 
support for the approximation ([m|) only in the case of ini- 
tially power-law spectra, we suggest that it is likely to be a 
good approximation for arbitrary power spectra. 



P!(k) = {21 + 1) / d s kP s {k)V e (^ k )W N {k,k) 



(67) 



The ratio Po(k)/P'2(k) of smoothed quadrupole-to- 
monopole harmonics satisfies the usual equation ( pCj ) in the 
linear limit. 

The whole hierarchy of smoothed spectra, for various 
N and £, can be derived from the generating function 



G(a,b) 



d 3 kP s (k)e 



-ak 2, — bk^fi^ 



(68) 



In terms of this generating function, the monopole harmonic 
of the redshift power spectrum, smoothed with the window 
W N is, 



Po(k) = 



■ d - 


N 


G(a,b = 0) 


. da_ 




_27rfc2iv+3 r [(2iV + 3)/2] 



(69) 



while the quadrupole redshift power smoothed with the same 
window is 



5 


■ d 


-3^~ 




■ d ~ 


AT-l 


2 




6 db_ 


6=0 


da_ 





G(a,b) 



27rfc2JV+3 r [(2AT + 3)/2] 



.(70) 



Substituting equation (^) for the redshift power spec- 
trum into the generating function 



G(a,b) = / d\A\e kq - k ^^ q) 



(71) 



3.5 The generating function of smoothed power 
spectrum multipoles 

So far we have obtained a number of analytic results for 
the redshift power spectrum in the Zel'dovich approxima- 
tion, but in general it is necessary to resort to numerics. 
Unfortunately, numerical integration of the redshift power 
spectrum, equation (|5^) , while feasible, is unpleasant (see 
Appendix A). One way to sidestep the numerical difficulties 
is to smooth the power spectrum. This is not such a bad idea 
because it is necessary to measure a smoothed power spec- 
trum from observations. That is, the Fourier modes mea- 
sured in a catalogue are already convolved with the window 
function of the catalogue, and one may choose to smooth 
the power spectrum further to reduce error bars. In addi- 
tion, smoothing makes defining quantities such a s th e zero 
crossing of the quadrupole, discussed in Section 3.6, more 
robust to random errors. 

The smoothed (unredshifted) power spectrum is defined 



P(k) 



d 3 kP(k)W N (k,~k) 



(65) 



We choose to adopt a power law times Gaussian smoothing 
kernel 



W N (k,k) 



fc 2Jv exp(-fc 2 /fc 2 ) 
27rfc2iv+3 r [(2iV + 3)/2)] 



(66) 



which is analytically convenient, and goes to a delta-function 
in the limit N — > oo. In subsequent Sections we use N — 1. 

For the redshift power spectrum, the smoothed harmon- 
ics can be defined similarly by 



(72) 



The k part of the integral ( |7l[ ) is the Fourier transform of a 
Gaussian, which solves in the usual way to 



f . , ^ 3/2 d 3 g 



9 2 *7, 1 /4 



(73) 



|$|l/2 

where |$| is the determinant of ^ij, 

|*| = (i> + + a){(^ + +i>-+a)[i> + (l+f3) 2 + a + b] 

+ /*V_[a(2/3 + (3 2 ) - b]} (74) 

with [i = z.q and tp±(q) given by equations ( pl| ) and (pj[), 
and = qi^^qj is the gg-component of the inverse of 

<- ^»*+« 

+ (l- f i 2 )ty + (2(3 + l3 2 )+b + v 2 i>-f3 2 ]} ■ (75) 

Evaluation of the smoothed monopole and quadrupole spec- 
tra from equations (^), (0) and (^) still involves a dou- 
ble integral over p and q. However, the integrand here is 
well-behaved (at least for small N, such as N — 1) and de- 
pends only on elementary functions (besides ij)±(q)). This 
compares to the unsmoothed case (Appendix A), where the 
integrand is ill-behaved, and evaluating the harmonics in- 
volves a double integral over an infinite sum involving spe- 
cial functions. 
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3.6 Numerical results for quadrupole-to-monopole 
power 

From the generating function of multipole moments, equa- 
tion ([73]), we have calculated numerically the smoothed 
quadrupole and monopole moments of the power spectrum 
for initially power law spectra Pl oc k n with various values 
of the index n, and for various f3. The smoothing kernel is 
Wjv, equation (pi]), with N = 1. Figure || shows the ratio R 
of the smoothed quadrupole-to-monopole power 



R : 



PS &) 



(76) 



for a representative sample of cases, n — —1, —1.5, —2, 
— 2.5 and —2.9, and for /3 = and 1 in each case. Values 
intermediate between these give curves intermediate to those 
plotted. The plotted ratio has been divided by the linear 
ratio, R L = |/3 2 )/(l + §/?+ \(3 2 ), equation @, and 

the wavenumber has been scaled to the zero-crossing scale, 
fc . 

It is immediately apparent from Figure ^| that the shape 
of the quadrupole-to-monopole ratio R depends on the spec- 
tral index n, but is insensitive to (3 at fixed n, at least for 
n i> —2.5. The insensitivity to f3 is an interesting result, and 
we have previously invoked this insensitivity in proposing 
the approximation ( |64| ) to the (unsmoothed) ratio. As to 
the spectral index, Figure ^ shows that for n ^ —2, the ra- 
tio R actually rises up above its linear value before turning 
over and going through zero. Physically, the linear squash- 
ing effect is first enhanced by caustics in redshift space as 
structures approach turnaround, and is then subsequently 
negated by "fingers-of-god" as structures collapse. Figure ti 
indicates that for n ^ — 2 the caustic effect wins over the 
finger-of-god effect at translinear scales. However, as the en- 
hancement peak scales as fc p k ~ 0.2 1 ^" +3 ' 1 ko, we find that 
R decreases for n Ss —2, for fixed /3, in the regime k ~ ko. 

A useful fitting function to the curves in Figure |^ is 

R = R L e~ {n+1)/2 [l - (fc/fc )°- 6( ™ +3) ] , (77) 

which is accurate to within 10 percent over the range 0.1 < 
k/ko < 1, and — 3 < n < — 1. We shall use the fitting func- 
tion ( fn| ) in fitting the observed ratio R measured in the 
QDOT plus 1.2 Jy survey to /3 and n, Section L2. 

The zero-crossing scale of the quadrupole power also 
contains useful information about (3 and the spectral index. 
The zero-crossing scale discussed below refers strictly to that 
of the smoothed quadrupole power, but the zero-crossing 
scale of the unsmoothed quadrupole is not greatly different. 

Suppose that the zero-crossing of the smoothed 
quadrupole power P| (&) is observed to occur at wavevector 
fco. Define a corresponding real scale by qo = a/ko, where 
a is a fitting constant, discussed below. Then we find that 
the dimensionless amplitude of the displacement correlation 
function at the zero-crossing scale go is fitted tolerably well 
by the following fitting formula 



1 + 13 



qo = a/ko 



(78) 



In the Zel'dovich approximation, formula (|78j) with the val- 
ues (we adopt here the value of a whi ch best fits the TV-body 
simulations reported in Section 4.1, since the same value 
works well also for Zeldovich) 




Figure 2. Plots of the ratio R of smoothed quadrupole-to- 
monopole power against (fc/fco) n + 3 for initially power-law spectra 
with indices, from lower curve to upper, of n = —1, —1.5, —2,-2.5 
and —2.9. In each case curves for (3 = and 1 are shown, which 
lie close to each other (and practically on top of each other for 
n = — 1 and —1.5) when R is scaled to the linear amplitude _Rl 
and the zero-crossing scale fco- For each value of the spectral in- 
dex, (3 = gives the slightly higher scaled ratio. 



a = 2.3, C* = 0.19 



(79) 



is accurate to better than 20% for initially power-law spec- 
tra over the observationally interesting range of spectral 
indexes — 2 < n < — 1 and linear growth rate parameter 
< /? < 1. As a matter of interest, in the case n = — 1 
the exact analytic result (in the Zeldovich approximation) 
is ■t/j+(<lo)/<lo = l/[o 2 (l + P)]- Physically, it is not too sur- 
prising that the zero-crossing of the quadrupole should oc- 
cur when the rms displacement 4'+ 2 (l) between pairs has 
reached an appreciable fraction of their separation q. 

The amplitude of the power spectrum itself at the zero- 
crossing fco is a more complicated function of spectral in- 
dex. During linear evolution, the dimensionless smoothed 
monopole redshift power Ao 2 (fc) is related to the dimen- 
sionless displacement correlation function ip+(<j)/<l 2 at q = 
a/k by (compare equation [Q; note the factor T[(5 + 
ri)/2]/F(5/2) below comes from smoothing) 



Ao 2 (fc) = {l + \l3+\p 2 )-{2-n)sm[{-n- 

6 5 7T 

r[(5+n)/2] a 3 +n V+(a/fc) 



l)7r/2]r(l-n) 



r(5/2) 



(a/fc) 2 



(80) 



Nonlinear evolution, which is starting to become important 
at zero-crossing, will modify the power from the linearly 
evolved value (BOh . Setting aside this effect of nonlinearity, 
one concludes from equation (|^), combined with the ap- 
proximation (fr8|), that the amplitude of the monopole power 
at zero-crossing is sensitive to the spectral index n, but rel- 
atively insensitive to f3 (the 1/(1 + f3) factor in equatio n |78| 
tending to cancel the 1 + 1/3 + \(3 2 factor in equation 



© 0000 RAS, MNRAS 000, 000-000 



10 A. N. Taylor & A. J. S. Hamilton 








Figure 3. Plots of the ratio R of smoothed monopole— to— 
quadrupole power scaled by the zero-crossing, ko. Prom top to 
bottom on scales k < ko are plotted the numerical results for 
the Zel'dovich evolution for each power spectrum (solid lines) for 
(3 = 1 with n = — 1, —1.5 and —2 and j3 = 0.5 with n = — 1, 
— 1.5 and —2. Overlayed are the results from a set of N-body sim- 
ulations with the same initial parameters (points) . Prrors on the 
simulations are from cosmic variance. 



It is worth emphasising this interesting if somewhat dis- 
appointing result: one might have hoped that the amplitude 
of the observed monopole power at the zero-crossing of the 
quadrupole would provide a measure of /3 in the mildly non- 
linear regime, offering a check on the value of j3 measured 
from the quadrupole-to-monopole ratio in the linear regime, 
and perhaps a test of the bias b. However, the hope fails: the 
power at zero-crossing is sensitive mainly to the spectral in- 
dex, not to (3. 



4 COMPARISON WITH SIMULATIONS AND 
OBSERVATIONS 

4.1 Comparison with simulations 

Given the approximate nature of the Zel'dovich mapping 
assumed in Sections |i| and it is important to check the 
accuracy of these results with TV-body simulations, where 
nonlinear effects can be accounted for on smaller scales. 
The simulations we ran used the AP 3 M code of Couchman 
(1991), generalised to low-f2 by Peacock & Dodds (1995), on 
a 64 3 grid with 64 3 particles. 

We ran simulations with power law initial spectra with 
indices n = — 1, —1.5 and —2, both for Q = 1, and for a 
low-f2 spatially fiat model with f2 = 0.3 and Q,a = 0.7, cor- 
responding to P = 0.5. In each case three simulations were 
computed with differing random Gaussian initial conditions, 
18 simulations altogether. We assumed no bias in the den- 
sity field, and evolved all the models for 300 time steps so 



as to erase the effects of using the Zel'dovich approximation 
in the initial setting up of the simulation. 

Power multipoles were calculated by assigning densities 
to the grid by the cloud-in-cell algorithm and FFT-ing the 
cube. The shot noise component of the power was subtracted 
and the resulting Fourier space compressed by averaging 
over azimuthal angles. The multipoles were then calculated 
by a least squares fit to each mode: 



A 2 (fc, Ml )-^A'(fc)^(^) 



(81) 



The weighting scheme used was the inverse variance, 
Wi(k) = l/a 2 (k, fii), calculated from the intrinsic cosmic 
variance, cr(fc) = \/2A 2 (fc), for each mode. The resulting 
power multipoles were then smoothed using the kernel func- 
tion Wn with TV = 1, equation do6|), and an inverse variance 
weighting scheme. 

Generally, we found that the Zel'dovich approximation 
ceased to provide a good approximation to the evolution 
of monopole or quadrupole power at moderately nonlinear 
epochs. This is to be expected, since as discussed in Section 



2.3 



the expansion of structures following collapse which oc- 
curs in the Zel'dovich approximation causes power to reach 
a maximum and then decline, whereas in reality power con- 
tinues to increase monotonically. 

However, as shown in Figure ^, the Zel'dovich approxi- 
mation does appear to reproduce remarkably well the ampli- 
tude and shape of the ratio R of the quadrupole-to-monopole 
power seen in the TV-body simulations, for scales down to the 
zero-crossing of the quadrupole, k < ko. Small deviations 
occur at small k, but these can be attributed to the lack 
of contribution from smaller wavenumbers in the smoothing 
process of the simulations. Small deviations are also seen in 
the n = — 2 spectra, at least some of which may arise from 
cosmic variance and the influence of the finite box size on 
the longest wavelengths in the simulations. But the change 
in shape of the quadrupole-monopole ratio with spectral 
index predicted by the Zel'dovich approximation is clearly 
seen in the simulations. All the models fail in the regime 
k > ko, where we expect virialised clusters to provide the 
strongest redshift distortions (Jackson 1972). 

The scaling relation suggested by equation j7q) for 
the amplitude of the displacement covariance at the zero- 
crossing of the quadrupole is also found to hold for the 
simulations. However, the amplitude C predicted by the 
Zel'dovich approximation was not so accurate. Applying a 
least squares fit of the model to the ensemble of 18 simula- 
tions, we found 



a = 2.3 ±0.4, C = 0.54 ±0.07 



(82) 



with a formal x = 15-6 for 16 degrees of freedom. In gen- 
eral while the form of the scaling relation (|78j) holds, the 
Zel'dovich approximation underestimates the wavenumber 
at zero-crossing of the quadrupole by a factor of approxi- 
mately 2 in the cases considered. This can be understood 
as arising from the outflow of particles which follows caus- 
tic formation in the Zel'dovich approximation. This out- 
flow exaggerates the "finger-of-god" effect, causing the zero- 
crossing of the quadrupole to appear at larger scales in the 
Zel'dovich approximation than is actually the case. 

In fitting to observed data in the next Section, we use 
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£/[h -1 Mpc] 



n 



Figure 4. Plot of the ratio R of smoothed monopole-to- 
quadrupole power as a function of wavenumber k for the merged 
QDOT plus 1.2 Jy redshift survey. The solid line is a fit to a 
model with f3 = 0.8 and n = —2.5, and with ko = 0.5/iMpc — . 



the iV-body simulations to calibrate the zero-crossing of the 
quadmpole. 



4.2 Application to the QDOT plus 1.2 Jy redshift 
survey 

We have applied the results obtained above to estimate the 
distortion parameter, /3, and the local primordial spectral 
index, n, from the amplitude and shape of the ratio R of 
smoothed quadrupole-to-monopole power observed in the 
merged QDOT plus 1.2 Jy IRAS redshift survey. We further 
infer the variance erf of counts in 8 ft~ 1 Mpc spheres from the 
observed scale of the zero-crossing of the quadrupole. The 
observed quadrupole-to-monopole ratio, taken from Hamil- 
ton (1995), is shown in Figure ^. Superimposed on the ob- 
servations is a model with f3 = 0.8 and n — —2.5. The 
zero-crossing is estimated to be at ko = 0.5 ± 0.1/iMpc -1 . 

We use the fitting function ((??]) to fit the amplitude, 
which depends on (3, and the shape, which depends on n, of 
the observed quadrupole-to-monopole ratio R. We construct 
a likelihood function 



L = exp —diMijiij 



(83) 



where = (data — model) is the difference between the 
i th data point and the model prediction, and M%j is the 
data covariance matrix whose diagonal is the uncertainty 
on each data point and whose off-diagonal terms are es- 
timated assuming the underlying data have independent 
Gaussian-distributed errors, and the major source of cor- 
relations comes from smoothing. As our model breaks down 
on scales smaller than the zero-crossing point, and the er- 
rors on the observed quadrupole diverge at low wavenum- 
ber, we restrict our analysis to wavemodes in the range 
0.06 < k/ihMpc- 1 ) < 0.5. 



Figure 5. Likelihood function C{j3, n|data), of fi and the local 
spectral index, n, from the quadrupole-to-monopole ratio of the 
merged QDOT plus 1.2 Jy redshift survey. Contours are spaced 
by SlnC = —0.5 with the inner contour delineating the 68% and 
outer contour the 95% confidence region. While the spectral index 
is not strongly constrained, /3 < 0.5 is ruled out at the 95% level. 
This rules out f2 < 0.3 if IRAS galaxies are unbiased. 



Figure g shows the marginalised likelihood in the (f3, n) 
plane. As n — > —3, j3 is forced to infinity by our fitting func- 
tion. The likelihood function selects this range as the better 
fit to the data due to the slight inflection in the quadrupole- 
monopole ratio around k — 0.2/iMpc" 1 . However if we con- 
sider only the 95% confidence limit given by the outer con- 
tour, low values of /3 < 0.5 are strongly excluded. In terms 
of the galaxy bias parameter, this implies that the density 
parameter Q > 0.36 5//3 . The value of the primordial spec- 
tral index on these scales is not strongly constrained, but is 
consistent with n ~ — 2. 

We estimate the clustering amplitude as from the scale 
of the zero-crossing of the quadrupole, using equation ( |T8| ) 
with the parameters (^) calibrated from the iV-body sim- 
ulations. In figure M we plot the marginalised likelihood 
function C(ag) = J~df3dnC(f3, n, ag|data). We deduce that 
erg = 0.7 ± 0.05 where we quote 95% errors. This is in very 
good agreement with values found by other methods for the 
IRAS galaxy catalogs (eg, Fisher et al. 1994, Heavens & Tay- 
lor 1995) if IRAS galaxies are unbiased. Given the indirect 
method of our measurement, we find this very encouraging. 



5 SUMMARY 

The main aim of this paper has been to study the redshift 
distortion of the power spectrum in the moderately nonlin- 
ear regime. Our approach is motivated by the consideration 
that structures may appear in redshift space more nonlinear 
than they really are — for example, regions which are just 
turning around in real space appear as caustics in redshift 
space. This leads us firstly to work in Lagrangian space, 
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Figure 6. Marginalised likelihood function C(ag) of the variance 
in 8 h~ 1 WLpc spheres inferred from the zero-crossing of the QDOT 
plus 1.2 Jy quadrupolc. 



and secondly to use the Zel'dovich approximation, which is 
in effect linear theory expressed in Lagrangian space. 

We started by deriving an expression relating the power 
spectrum to the Lagrangian displacement field. We used 
this to determine the evolution of the power spectrum in 
the Zel'dovich approximation, first in real space, then in 
redshift space in the plane-parallel approximation. We pre- 
sented some analytic solutions for initially power-law spec- 
tra. In particular, we showed that a spectrum with index 
n = — 1 gives rise to an incoherent Gaussian displacement 
field, producing nonlinear redshift distortions with the same 
form as the Kaiser-Peacock model. 

We derived various analytic and numerical results for 
the observationally interesting ratio R of quadrupole-to- 
monopole redshift power, whose value in the linear regime, 
#l = (§/3+f/3 2 )/(l + §/3+±/3 2 ), yields ameasure of the dis- 
tortion parameter f3. In the Zel'dovich approximation, the 
amplitude of the ratio R is set by the linear value J?l, but the 
shape of J? as a function of wavenumber k depends mainly on 
the spectral index n, and is insensitive to /3 at fixed n. The 
zero-crossing of the quadrupole power occurs at the point 
where the Zel'dovich power spectrum is a maximum as a 
function of time. 

We have tested the Zel'dovich results against iV-body 
simulations with initially power law spectra, in both Q — 1 
and low-fi models. The simulations show that the Zel'dovich 
approximation ceases to provide a good approximation to 
the power spectrum at moderately nonlinear epochs. Re- 
markably however, the Zel'dovich approximation predicts 
rather well the amplitude and shape of the quadrupole-to- 
monopole ratio R on scales down to the zero-crossing of 
the quadrupole, when R is scaled to the scale of the zero- 
crossing. The Zel'dovich approximation underestimates the 
wavenumber at zero-crossing by a factor of about two, al- 



though it predicts correctly the way the zero-crossing scales 
with (3 and n. 

We have applied these findings to estimate the distor- 
tion parameter /3, the local spectral index n, and the vari- 
ance a\ of counts in 8 /i~ 1 Mpc spheres, from the quadrupole- 
to-monopole ratio R measured in the merged QDOT plus 1.2 
Jy redshift survey. We find that the distortion parameter is 
constrained to /3 > 0.5 at the 95% level. The spectral in- 
dex is not well constrained, but is consistent with n £s — 2 
at translinear scales. The clustering amplitude as, inferred 
from the scale of the zero-crossing of R, is ag = 0.7 ± 0.05, 
consistent with other estimates. 

The success of the Zel'dovich model in describing the 
quadrupole-to-monopole ratio R suggests that departures 
from the linear value _Rl at translinear scales are caused 
mainly by infall on to clusters, not by virialised cluster cores. 
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APPENDIX A 

This Section shows how to evaluate the 'difficult' integral 
( [53| ) for the redshift power spectrum in the Zel'dovich ap- 
proximation. The integral is a two-dimensional integral over 
H — K.q and separation q. We choose to develop the inte- 
gral over [s, analytically as an infinite sum, both because the 
numerical integral over jj, is itself unpleasant, and because 
the expression as an infinite sum allows the q-integrand to 
be continued analytically into the complex g-plane, which 
stabilizes the integration over q. 



y (ATV-)" ( d \ 2 



For small kqs, the Bessel function Jo{kqs*Jl — /J?) in 
the integrand can be expanded as a power series in 1 — fi 2 . 
Integration of the leading term (which is 1) in the series 
leads to an expression similar to and higher terms in 

1 — fi 2 can be generated by repeated differentiation of this 
expression. The result is the infinite sum 



^ n! \dAj A=kqc i[Ai + (kqs)2]i/2 

n—Q 



(85) 



Expression ( pq ) works when K 2 ip^ is small and kqs is not 
small, which is complementary to the domain of validity of 
equation (|S4|). A program such as Mathematica helps con- 
siderably with the algebra involved in taking the repeated 
derivatives in both equations (|84h and (B5I). 



The integration of equations ( |84j ) and ( |85[ ) over q must 
be done numerically. As in the unredshifted case, the oscilla- 
tory character of the integrand for large kq can be converted 
into exponential decay by shifting the path of integration of 
q into the upper complex plane. One should be careful here 
that although the real parts of the integrands of equations 
( [si] ) and ( [gq ) agree for real q, their imaginary parts differ, so 
that their analytic continuations in the complex q-plane dif- 
fer. The integrands could be made identical by adding their 
complex conjugate expressions, but this would introduce an 
undesirable exponentially increasing component in the up- 
per complex plane. Therefore, if it is necessary to use both 
expressions (g^) and (0) in integrating from q = to oo, 
then the transition from one expression to the other along 
the path of integration must occur on the real line. 

In the cases we tried, an integration path similar to 
that described at the end of Section 2.2 proved satisfac- 



tory. Namely, integrate first some way along the real axis, 
then complete the integration along a straight line upwardly 
slanted in the complex plane. In the case K 2 ip-(q = 1/ks) > 
1, equation ( [si] ) works everywhere on the path of integra- 
tion. In the case K 2 ip-(q — 1/ks) < 1, equation (^J) works 
for the first part of the integration, along the real line, but 
equation (pa) is required for the second part, a suitable path 
being vertically upwards in the complex plane. 

Computing the quadrupole and monopole harmonics of 
the redshift power spectrum requires a final numerical in- 
tegration of P s {k) over fik with Vt{^h) for I = and 2, 
equation (Ba). 



P s (k) =TZe I 2Tvq 2 dqe 



E 



(-r(kqs/2) 2n / d y 

n! 2 \dBj B=K 2. 



F(kqc,B) (84) 



where F is the function defined by equation (|26|) , and s and 
c are given by equations ( |54| ) and (^H|) . The sum in equation 
( p4| ) can be truncated to a few terms if kqs is small. For large 
k (nonlinear limit), this truncation effectively works for all 
q, since then the exponential factor e~ K causes the 

integrand to become negligible when kqs is not small. 



The expansion (|34|) fails when kqs is large and at the 
same time K 2 + if)_ ) is small. Here a series expansion in 
K 2 ip- leads to the following infinite sum: 



P 8 (k) =TZe 4Tvq 2 dqe 
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